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Abstract: We revisit the fundamental topic of light scattering by single 

homogenous nanoparticles from the new perspective of excitation and ma¬ 
nipulation of toroidal dipoles. It is revealed that besides within all-dielectric 
particles, toroidal dipoles can also be efficiently excited within homogenous 
metallic nanoparticles. Moreover, we show that those toroidal dipoles 
excited can be spectrally tuned through adjusting the radial anisotropy 
parameters of the materials, which paves the way for further more fiexible 
manipulations of the toroidal responses within photonic systems. The study 
into toroidal multipole excitation and tuning within nanoparticles deepens 
our understanding of the seminal problem of light scattering, and may 
incubate many scattering related fundamental researches and applications. 

© 2015 Optical Society of America 

OCXS codes: (290.5850) Scattering, particles; (290.4020) Mie theory; (160.1190) Anisotropic 
optical materials 


References and links 

1. T. Kaelberer, V. Fedotov, N. Papasimakis, D. Tsai, and N. Zheludev, “Toroidal dipolar response in a metamate¬ 
rial,” Science 330 , 1510 (2010). 

2. Z. Dong, P. Ni, J. Zhu, X. Yin, and X. Zhang, “Toroidal dipole response in a multifold double-ring metamaterial,” 
Opt. Express 20 , 13065 (2012). 

3. Y.-W. Huang, W. T. Chen, P. C. Wu, V. Fedotov, V. Savinov, Y. Z. Ho, Y.-F. Chau, N. I. Zheludev, and D. P. Tsai, 
“Design of plasmonic toroidal metamaterials at optical frequencies,” Opt. Express 20 , 1760 (2012). 

4. B. Ogut, N. Talebi, R. Vogelgesang, W. Sigle, and P. A. van Aken, “Toroidal plasmonic eigenmodes in oligomer 
nanocavities for the visible,” Nano Lett. 12 , 5239 (2012). 

5. V. A. Fedotov, A. Rogacheva, V. Savinov, D. Tsai, and N. I. Zheludev, “Resonant transparency and non-trivial 
non-radiating excitations in toroidal metamaterials,” Sci. Rep. 3, 2967 (2013). 

6. Y. Fan, Z. Wei, H. Li, H. Chen, and C. M. Soukoulis, “Low-loss and high-q planar metamaterial with toroidal 
moment,” Phys. Rev. B 87, 115417 (2013). 

7. V. Savinov, V. Fedotov, and N. Zheludev, “Toroidal dipolar excitation and macroscopic electromagnetic proper¬ 
ties of metamaterials,” Phys. Rev. B 89, 205112 (2014). 

8. Q. Zhang, J. J. Xiao, X. M. Zhang, D. Han, and L. Gao, “Core-shell-structured dielectric-metal circular nanodisk 
antenna: Gap plasmon assisted magnetic toroid-like cavity modes,” ACS Photon. 2, 60 (2014). 

9. A. A. Basharin, M. Kafesaki, E. N. Economou, C. M. Soukoulis, V. A. Fedotov, V. Savinov, and N. I. Zheludev, 
“Dielectric metamaterials with toroidal dipolar response,” Phys. Rev. X 5 , 011036 (2015). 

10. A. E. Miroshnichenko, A. B. Evlyukhin, Y. E. Yu, R. M. Bakker, A. Chipouline, A. I. Kuznetsov, B. Lukyanchuk, 
B. N. Chichkov, and Y. S. Kivshar, “Seeing the unseen: observation of an anapole with dielectric nanoparticles,” 
Nat. Commun. (2015). arXiv: 1412.0299 



11. W. Liu, J. Zhang, and A. E. Miroshnichenko, “Toroidal dipole induced transparency for core-shell nanoparticles,” 
Laser Photon. Rev.,doi:10.1002/lpor.201500102 (2015). arXiv: 1412.4931 

12. R. M. Bakker, D. Permyakov, Y. F. Yu, D. Markovich, R. Paniagua-Domlnguez, L. Gonzaga, A. Samusev, 
Y. Kivshar, B. Lukyanchuk, and A. I. Kuznetsov, “Magnetic and electric hotspots with silicon nanodimers,” 
Nano Lett. 15 , 2137 (2015). 

13. W. Liu, J. Zhang, B. Lei, H. Hu, and A. E. Miroshnichenko, “Invisible nanowires with interfering electric and 
toroidal dipoles,” Opt. Lett. 40 , 2293 (2015). 

14. V. M. Dubovik and V. V. Tugushev, “Toroidmoments in electrodynamics and solid-state physics,” Phys. Rep. 
187 , 145 (1990). 

15. E. E. Radescu and G. Vaman, “Exact calculation of the angular momentum loss, recoil force, and radiation 
intensity for an arbitrary source in terms of electric, magnetic, and toroid multipoles,” Phys. Rev. E 65 , 046609 
( 2002 ). 

16. J. D. Jackson, Classical electrodynamics (Wiley New York, 1962). 

17. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983). 

18. L. Novotny and B. Hecht, Principles of nano-optics (Cambridge University Press, Cambridge, 2012), 2nd ed. 

19. W. Liu, R. F. Oulton, and Y. S. Kivshar, “Geometric interpretations for resonances of plasmonic nanoparticles,” 
Sci. Rep. 5 , 12148 (2015). 

20. A. I. Kuznetsov, A. E. Miroshnichenko, Y. H. Fu, J. B. Zhang, and B. S. Lukyanchuk, “Magnetic light,” Sci. Rep. 
2 , 492 (2012). 

21. A. B. Evlyukhin, S. M. Novikov, U. Zywietz, R. L. Eriksen, C. Reinhardt, S. I. Bozhevolnyi, and B. N. Chichkov, 
“Demonstration of magnetic dipole resonances of dielectric nanospheres in the visible region,” Nano Lett. 12 , 
3749 (2012). 

22. W. Liu, A. E. Miroshnichenko, and Y. S. Kivshar, “Control of light scattering by nanoparticles with optically- 
induced magnetic responses,” Chin. Phys. B 23 , 047806 (2014). 

23. W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B 39 , 9852 (1989). 

24. M. S. Wheeler, J. S. Aitchison, and M. Mojahedi, “Coated nonmagnetic spheres with a negative index of refrac¬ 
tion at infrared frequencies,” Phys. Rev. B 73, 045105 (2006). 

25. I. Fernandez-Corbaton, S. Nanz, and C. Rockstuhl, “On the dynamic toroidal multipoles,” arXiv: 1507.00755 
(2015). 

26. I. Fernandez-Corbaton, S. Nanz, R. Alaee, and C. Rockstuhl, “Exact dipolar moments of a localized electric 
current distribution,” arXiv: 1507.00752 (2015). 

27. C. C. Lam, P. T. Leung, and K. Young, “Explicit asymptotic formulas for the positions, widths, and strengths of 
resonances in mie scattering,” J. Opt. Soc. Am. B 9 , 1585 (1992). 

28. W. Liu, “Ultra-directional super-scattering of homogenous spherical particles with radial anisotropy,” Opt. Ex¬ 
press 23 , 14734 (2015). 

29. P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370 (1972). 

30. C.-W. Qiu, L.-W. Li, T.-S. Yeo, and S. Zouhdi, “Scattering by rotationally symmetric anisotropic spheres: Poten¬ 
tial formulation and parametric studies,” Phys. Rev. E 75 , 026609 (2007). 

31. C. Qiu, L. Gao, J. D. Joannopoulos, and M. Soljacic, “Light scattering from anisotropic particles: propagation, 
localization, and nonlinearity,” Laser Photonics Rev. 4 , 268 (2010). 

32. Y. X. Ni, L. Gao, A. E. Miroshnichenko, and C. W. Qiu, “Controlling light scattering and polarization by spherical 
particles with radial anisotropy,” Opt. Express 21 , 8091 (2013). 

33. W. Liu, A. E. Miroshnichenko, D. N. Neshev, and Y. S. Kivshar, “Broadband unidirectional scattering by 
magneto-electric core-shell nanoparticles,” ACS Nano 6, 5489 (2012). 

34. W. Liu, J. Zhang, B. Lei, H. Ma, W. Xie, and H. Hu, “Ultra-directional forward scattering by individual core-shell 
nanoparticles,” Opt. Express 22 , 16178 (2014). 

35. R. F. Oulton, V. J. Sorger, T. Zentgraf, R. M. Ma, C. Gladden, L. Dai, G. Bartal, and X. Zhang, “Plasmon lasers 
at deep sub wavelength scale,” Nature 461 , 629 (2009). 

36. M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Sutee- 
wong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460 , 1110 (2009). 

37. A. V. Kabashin, P. Evans, S. Pastkovsky, W. Hendren, G. A. Wurtz, R. Atkinson, R. Pollard, V. A. Podolskiy, and 
A. V. Zayats, “Plasmonic nanorod metamaterials for biosensing,” Nat. Mater. 8, 867 (2009). 

38. L. Novotny and N. van Hulst, “Antennas for light,” Nat. Photon. 5, 83 (2011). 

39. H. A. Atwater and A. Polman, “Plasmonics for improved photovoltaic devices,” Nat. Mater. 9, 865 (2010). 

40. L. Lu, J. D. Joannopoulos, and M. Soljacic, “Topological photonics,” Nat. Photon. 8, 821 (2014). 

41. F. Xia, H. Wang, D. Xiao, M. Dubey, and A. Ramasubramaniam, “Two-dimensional material nanophotonics,” 
Nat. Photon. 8, 899 (2014). 

42. G. N. Afanasiev and Y. P. Stepanovsky, “The electromagnetic field of elementary time-dependent toroidal 
sources,” J. Phys. A Math. Gen. 28, 4565 (1995). 

43. E. A. Marengo and R. W. Ziolkowski, “Nonradiating sources, the aharonov-bohm effect, and the question of 
measurability of electromagnetic potentials,” Radio Sci. 37, 19 (2002). 



1. Introduction 


Stimulated by the recent experimental demonstration of toroidal dipoles (TDs) excited on the 
platform of metamaterials consisting of specially organized split ring resonators d, studies into 
TDs and more generally toroidal multipoles (TMs) within photonic nano-systems have attached 
surging attention lEO . Compared to conventional electric and magnetic multipoles, TMs can 
be viewed as poloidal currents flowing on the surface of a torus along its meridians, or as a series 
of magnetic dipoles (enclosed circulating currents) aligned along enclosed paths IMKTSI. For 
the near held in terms of charge-current distributions, TMs are contrastingly different from their 
electric and magnetic counterparts Eni. In the far held however, TMs and conventional electric 
multipoles are indistinguishable Bin, which results in the situation that in various scattering 
theories concerning particle scattering, TMs are not considered as separate entities ifT^fTTl . 
Generally under most circumstances, the contributions of TMs have been incorrectly attributed 
to conventional electric multipoles, rendering TMs unfortunately as missing or hidden terms. 

Following the first experimental demonstration of TDs (H, TMs have been observed in var¬ 
ious configurations, most of which (influenced by the configuration designed in Ref. ID) are 
kind of toroid-shaped composite structures ||2]-[6l[9]l. This is understandable: the corresponding 
charge-current distribution of TMs is in a toroidal shape and naturally in such composite struc¬ 
tures TMs can be efficiently excited. But this unfortunately leads to more or less a dogma that 
TMs can only be observed in carefully engineered compound structures and consequently there 
have been very few searches for the existence of TMs in simple individual structures il. Only 
recently it is shown decisively that TMs can be excited within single scattering particles of 
fundamental (spherical or cylindrical) shapes, which also have been proven to play indispens¬ 
able roles for some exotic scattering properties |[T0l[TTl[T3. The incorporation of TMs extends 
significantly our understanding of the fundamental problem of light scattering by spherical and 
cylindrical particles. Moreover, considering that the scattering of spherical and cylindrical par¬ 
ticles is the cornerstone of many scattering related fundamental studies and applications, the 
incorporation of TMs into scattering photonic systems can possibly spur further breakthroughs 
in sensing, nanoantennas, nanoscale lasing and photovoltaic devices ifT^ . 

Here in this article we restudy the seminal problem of plane wave light scattering by homoge¬ 
nous spherical nanoparticles, but from a new angle of excitation and manipulations of TDs. It 
is demonstrated that, similar to all-dielectric particles, even a homogenous metallic spherical 
nanoparticle can surprisingly support TDs in the spectral region beyond the quasi-static regime. 
We further show that radial anisotropy can be employed to tune spectrally the TDs excited, pro¬ 
viding extra freedom for TMs manipulations, which can possibly And significant applications 
in TMs related topics, especially those where near-held distributions of fields and currents play 
a crucial role. 

2. Methods and Expressions 

2.7. Conventional Multipole Expansions and Emerging Toroidal Multipoles 

To obtain the radiated fields of an arbitrarily distributed charge-current source, the general 
approach is to firstly calculate the retarded vector potentials and then get the associated electro¬ 
magnetic fields from those potentials . When the retardation effect within the charge-current 
source can be neglected (in other words the geometrical dimension of the source is far smaller 
than the effective wavelength and thus the quasi-static approximation can be applied), the vec¬ 
tor potentials can be Taylor-expanded into multiple terms of different orders, from which the 
conventional electric and magnetic multipoles can be deduced M- For example, the zeroth 
order expansion term corresponds to the conventional electric dipole P [the exp(—icot + ik- r) 
notation has been adopted where k and co denote wave-vector and angular frequency respec- 



Fig. 1. ((a) Schematic illustration of the scattering of an incident plane wave by a spherical 
particle of radius R. The spherical particle can be homogenous and isotropic with reflective 
index n, or can be radially anisotropic with radial permittivity £r and transverse permittivity 
St =n^- The anisotropy parameter is defined as T] = (r/ = 1 for isotropic spheres). The 

plane wave is polarized (in terms of electric field) along x direction and is propagating along 
z direction, (b) Schematic illustration of the toroidal dipole excitation within the spherical 
particle. Both the current J and magnetic field H distributions have been shown. 


tively]: 

P=-L fd^rJir), (1) 

— ICO J 

and the first order expansion term includes both conventional electric quadrupole and magnetic 
dipole M: 

M=T/j3^[r xj(r)], (2) 

where c is the speed of light. Since M corresponds to the higher order term of the Taylor 
expansion, it results in more or less the mind-sets that: (1) P should be stronger than the M and 
(2) like most higher order modes O, M should be supported at higher frequencies than P. 
Nevertheless, we should keep in mind that the Taylor-expansion of the vector potentials ifT^ is 
based on the quasi-static approximation, and thus when the precondition of this approximation 
does not exist, the above conclusion is not valid anymore. An outstanding example of this is the 
recent demonstration of optically-induced magnetic response of higher permittivity dielectric 
particles ll20H22ll . where M is supported at lower frequencies and can be stronger than the 
electric dipole. 

The other consequence of the breakdown of the quasi-static approximation is that the con¬ 
ventional electric and magnetic dipoles [Eqs. ([T])-@] and higher order multipoles do not form 
a complete set and thus other multipoles have to be included as correcting terms ^T WTaWTSI . 
For dipolar responses, besides M and P, the lowest order correcting term would be the toroidal 
dipole: 

T=^/fifV[(r-J(r))r-2r2j]. (3) 

The scattered power of P and T are respectively: 


Wp 


^|p|2 

line ' ' 


Wt 


|2 

line ' ' 


( 4 ) 


where /Tq is the vacuum permeability. 










2.2. Electric and magnetic multipoles associated with spherical harmonics of Mie scattering 
particles 

For spherical (both homogenous and multi-layered) particles with incident plane waves of elec¬ 
tric field Eq, the scattered electric field can be expressed as ifTTll 


Es(r) - Eyyi\iayf^T^{v) 5 

m=l 


(5) 


where Um and bm are the Mie scattering coefficients; and Nm are vector spherical harmon- 
ics;and£;„ = /'”£o;;^ . In the far field the scattered fields [Eq.Q for r ^ oo] associated with 
am {Emiam^m{^)] and hm {EmhmM^m{^)] are exactly the same as the far fields of conventional 
electric and magnetic multipoles of order m, with conventional electric and magnetic dipoles 
P and M (m = 1) defined in Eqs.([T])-(|2l). For this far-field equivalence, it is natural to define 
another spherical harmonics based electric dipole as |[IOl[IIl[nil23ll2l: 


P((2i) = £o—p—Eo, 


( 6 ) 


where £o is the vacuum permittivity, and the associated scattering power of P(ai) is HIQIII111151 : 

|2 


Wp(aO = 


3k\Eo\^ 

kcopo ’ 


(7) 


We note here that F{ai) has been obtained based on the far scattered field expansion and thus 
has included all the dipolar transverse magnetic components HU El- In contrast, the con¬ 
ventional electric dipole P shown in Eq.([T]) has been obtained based on the near-field current 
integration, and here the current J(r) can be expressed as: 


J(r) = -/meo[e(r)-l]E(r), (8) 

where e(r) is a scalar if the particle is isotropic and a tensor if the particle is radially anisotropic. 
The conventional electric dipole P is only the lowest order expansion term of P(ai) shown in 
Eq.®. Similarly, the toroidal dipole T shown in Eq.® is the next higher order expansion of 
P(ai). So basically V{ai) can be expressed as ifTOlfTTlfTHfTSl : 

p(ai)=P + /^T+it, (9) 

where it denote higher order correcting terms of P(ai) with respect to kR Ill4[|15[[25l[26ll . As is 
shown in Eq.® the multipole expansion of the far scattered fields based on spherical harmonics 
is complete [Eqs.®-®] while the conventional multipole expansion based on the integration 
of the near-field current is incomplete [Eqs. O-©]. As the normalized particle radius (a = kR) 
is getting larger and larger (that is to say the particle radius R is getting larger and larger or the 
functioning wavelength is getting smaller or smaller), higher order correcting terms including 
the toroidal dipole T would arise. That is to say, multipole expansions even including toroidal 
multipoles are still incomplete. For sufficiently larger a, other higher order correcting terms 
must also be incorporated ifTOllllillSUTSlI^I^ . 


3. Results and Discussions 

In Fig. (Ha) the schematic of the scattering configuration is shown: the pane wave is scattered 
by a spherical particle, which could be homogenous and isotropic (metal or dielectric) with re¬ 
fractive index n or radially anisotropic with radial permittivity £r and transverse (perpendicular 










Fig. 2. (a) Scattered power spectra for a dielectric {n = 3) sphere (inset) of R = 200 nm. 
The contributions from current-integrated [based on Eq. 0 and Eq.©] electric dipole (P, 
red curves) and TD (T, blue curves), and those from the far-held deduced dipole [P(ai), 
dashed black curves] are shown. Eour points of different a = kR are indicated by black 
dots in (a) (A: a = 0.8, B: a = 1.5, C: a = 2.6, D: a = 3). The the corresponding near 
helds on the v — z plane of y = 0 are shown respectively in (A)-(D) in the bottom row. 
The distributions for both the normalized magnetic held along y direction (color-plots) 
and displacement held D = eE on plane (vector-plots, only the held inside the particle is 
plotted) are shown (the dashed black lines denote the boundaries of the particles), as is also 
the case for Pigs. [3P1 


to the radial direction) permittivity £f = n^. The anisotropy parameter is dehned as 7] —E^jEy. 
The plane wave is polarized along x direction and propagating along z direction. In Fig. \V[h) 
we show the excited toroidal dipole T within the particle by the magnetic held H and current 
distribution J. 

3.1. TD excitation within homogenous and isotropic dielectric particles 

As a hrst step we study the plane wave scattering by a homogenous and isotropic dielectric 
sphere of R = 200 nm and n = 3. The scattered power spectra [in terms of P(ai), P and T] are 
shown in Fig. ©a). Four points A-D are selected and indicated by black dots in Fig. ©a) and the 
corresponding near-held distributions are shown in Fig. ©A-D): the vector-plots correspond to 
the displacement held D = eE on the x — z plane [only the held within the particle is plotted, as 

















































Fig. 3. (a) Scattered power spectra for a silver sphere (inset) of R = 200 nm. The contribu¬ 
tions from P, T, and P(ai) are shown. Two points of different a are indicated by black dots 
in (a) (A: a = 0.8 and B: a = 3.5). The the corresponding near fields are shown in (A)-(B) 
on the right. 


is also the case for Fig. Q-Fig. (O] and the color-plots correspond to the normalized magnetic 
fields along y direction (H^). 

At point A of a = 0.8, the quasi-static approximation can be applied, and thus the total 
dipolar field can be approximated as a conventional electric dipole: P(< 2 i) ~ P. This is due 
to the fact that in the quasi-static spectral regime, T and other higher order correcting terms in 
Eq. Q can be neglected. To further confirm this, we show the near fields at point A in Fig. [2 A), 
where a typical electric dipole field is exhibited. 

At Points B and C of a = 0.8 and 2.6 respectively, which is beyond the quasi-static spectral 
regime, the quasi-static approximation breaks down and thus the contributions of TD have to 
be included [see Eq. ®]. At those points, P been totally suppressed and there are only effi¬ 
cient TD excitation: P(ai) ^ ikT. The corresponding near-field distributions are shown respec¬ 
tively in Fig.[2B-C). It is worth mentioning that the field distribution at points B in Fig.[2B) 
corresponds exactly to the filed and currents shown in Fig. [Jb), which is a typical and fun¬ 
damental representation of TD [the magnetic field is confined within the particle and there 
are two field intensity peaks (thus two oppositely circulating current loops on plane) across]. 
The field distribution at point C shown in Fig. [2C), though different from what is shown in 
Fig.EB), also corresponds to a TD, which however is of a higher mode number with more 
field lobes |[T9l[27l[28l [the magnetic field is also confined within the particle while there are 
four field intensity peaks across]. 

At Points D of a = 3 which is far beyond the quasi-static spectral regime, even the incorpo¬ 
ration of TD will not account for all the dipolar scattering, which can be confirmed in Fig. [2a): 
though at point D the overall dipolar scattering is negligible [P{ai), dashed black curve], TD 
response T is dominant over P and this means that they cannot interfere totally in a destruc¬ 
tive way, thus catcalling the scattering of each other HI]US; As a result, besides P and T, 
there must be higher order correcting terms for P{ai) as shown in Eq. (|9]). This further proves 
that the multipole expansions based on near-field integration is still incomplete even of the 
toroidal responses have been considered. Besides conventional electric, magnetic and toroidal 
multipoles, there are other terms which will become more and more significant as the particle 
dimension gets larger and larger compared to the functioning wavelength ([TdldSlEIlEll. The 
corresponding near-field distribution of point D is shown in Fig.|2T)). 































Fig. 4. Scattered power spectra for dielectric spheres (inset) of radius R = 200 nm with 
anisotropy parameters r] = £t/£r = I/3 in (sl) and T] = 2 in (b). The transverse permittivity is 
fixed at £t = 9. The contributions from P, T, and P(ai) are shown. Four points are indicated 
by black dots in (a) and (b) (A: a = 1.27, B: a = 2.31, C: a = 1.72, D: a = 2.8). The 
corresponding near fields are shown respectively in (A)-(D) on the right. 


3.2. TD excitation within homogenous and isotropic silver particles 

As a next step we study the plane wave scattering by a homogenous and isotropic Ag sphere of 
R = 200 nm and the permittivity of silver is adopted from the experimental data in Ref. Il29ll . 
The scattered power spectra [in terms of P(ai), P and T] are shown in Fig.^a). Two points 
A and B are selected and indicated by black dots in Fig. [Sja) and the corresponding near¬ 
field distributions are shown in Fig. OA-B). Similar to the dielectric sphere discussed in the 
subsection above: At point A of a = 0.8 in the quasi-static regime, P is the only present term 
[P(( 2 i) ^ P, see Fig.Oa)] and the near field exhibits a conventional electric dipole configuration 
[see Fig.[3tA)]; At point B of a = 3.5 beyond the quasi-static spectral regime, the TD term 
arises [see Fig. Oa)] and the corresponding near field shown in Fig. OB) is similar to that 
shown in Fig. 12B). The difference is that in Fig. [2B) there is only TD excitation while in 
Fig. OB) there are also P and higher order correcting terms, similar to point D in Fig. [2D). 
In a word, here we show that for homogenous metallic particles, T can be effectively excited 
beyond the quasi-static spectral regime, which can be stronger over P [Fig. [3a)]. 

3.3. Tuning the TDs excited within homogenous dielectric particle through radial anisotropy 

At the end, we study the scattering of radially anisotropic dielectric spheres of radius R = 
200 nm. The anisotropy parameter is 7] = £?/Cr and the transverse permittivity is fixed at = 9. 














































For radially anisotropic particles, both the Mie scattering coefficients and the near-fields can 
be calculated analytically in a similar manner to that of isotropic particles through modifying 
the expressions of the orders of related spherical Bessel functions to include the anisotropy 
parameters ll^[30H32ll . As a result, current J [as is shown in Eq. ^ where e(r) now is a 
tensor], P, T, P(ai) and their scattered power for radially anisotropic spherical particles can 
also be calculated directly. 

The scattered power spectra are shown in Fig. HJa) of r/ = 1/3 and in Fig. (Hb) of r/ = 2 
with contributions from P, T, and P(ai). Four points where only TD are excited [P(< 2 i) ~ 
ikT] with both suppressed P and higher order correcting terms are indicated by black dots (A: 
a = 1.27, B: a = 2.31, C: a = 1.72, D: a = 2.8). The corresponding near fields are shown 
respectively in Fig.[4tA)-(D). It is known that through employing radial anisotropy, P(ai) can 
be tuned to larger/smaller wavelengths (smaller/larger a) through employing smaller/larger 
rj (281. Similarly, TD excited within dielectric spheres can also be tuned by adjusting the radial 
anisotropy parameters, which is clear according to Fig.^a) and (b). Compared to the isotropic 
case of 7] = 1 shown in Fig. [2a), the TD responses [e.g. characterized by the positions where 
only TD is excited, as indicated by points B and C in Fig.|2a) and points A-D in Fig.[2a-b)] can 
be tuned to larger wavelengths (smaller a) through making 7] < 1 [see Fig.lUa)] and to smaller 
wavelengths (larger a) through rendering 77 > 1 [Fig. [dtb)]. For the near-field distributions, 
both TD with low mode number [see Fig. [2 A) and (C), which is similar to Fig. |2B)] and 
higher mode number [see Fig. (UB) and (D), which is similar to Fig. |2C)] can be effectively 
tuned within radially anisotropic dielectric particles. 

4. Conclusions and Outlook 

To conclude, we revisit the fundamental problem of Mie scattering of spherical particles from 
the new perspective of TD excitation and tuning. Besides the demonstration of TDs within 
dielectric particles, we further demonstrate that TDs can be effectively excited even within 
homogenous metallic particles as long as it is beyond the quasi-static spectral regime. It is 
further shown that the TDs excited can be spectrally tuned by adjusting the radial anisotropy 
parameters, which paves the way for further fiexible manipulations of toroidal responses. 

It is worth mentioning that here we have confined our study to the electric dipolar responses 
(transverse magnetic dipolar responses (111I17II ) and similar studies can certainly conducted for 
the magnetic dipolar responses (transverse electric dipolar responses GHEl), or both trans¬ 
verse magnetic and transverse electric responses of higher orders. Such multipole expansions 
including toroidal responses can also be implemented for other non-spherical structures where 
the fields and currents, and thus related multipoles can be calculated through numerical meth¬ 
ods. As for the expansion of the dipolar scattering shown in Eq. O, we have only discussed in 
detail the correcting terms up to toroidal dipoles, and similarly higher order correcting terms 
can also be further investigated |[TQl[TTl[T3tiT5l[^[^ . For the tuning of TDs excited. Here 
in this paper we employ only the radial anisotropy, and it is natural to expect that the other 
kinds of anisotropy, nonlinear and/or gain effects can also be used for further more fiexible 
manipulations of toroidal multipole responses. 

We believe that the incorporation of toroidal multipoles into scattering nanoparticles and the 
explorations into their tunability have not only expanded our understanding of the fundamen¬ 
tal Mie scattering physics, but also would be of great significance for both far-field scattering 
pattern shaping (22l|33l|34]| and near-field manipulations, which could possibly stimulate many 
scattering related fundamental researches and applications such as nano-lasing ||35l|36||, sens¬ 
ing [371, nanoantennas (381, and photovoltaic devices (39l . Moreover, investigations into light- 
matter interactions involving toroidal multipoles can be possibly conducted on the recently 
emerging topological (4Ql and/or two-dimensional (411 photonic platforms, and probably can 




lay a solid foundation for the decisive solution of the dynamical Aharonov-Bohm effect con¬ 
troversy concerning quantum information processing II421I43II . 
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